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A pointlike electric charge close to the surface of a three dimensional topological insulator induces 
a magnetic monopole mirror charge. We study the distance dependence of the statistical angle 
describing this induced dyon system. We find that the total angular momentum, which sometimes 
is used as signature of the statistical angle, for an electron outside a finite size spherical or tube 
shaped topological insulator is zero for any finite distance between the electron and the surface. 
However, we show that in the 2-electron system one can indeed isolate a non-trivial statistical 
angle for intermediate size loops, that is loops much larger than the distance of the charge to the 
sample but much smaller than the size of the sample. The necessity for this limit confirms the 2+1 
dimensional nature of the non-trivial statistical angle found in previous work. Our results clarify the 
conditions under which the statistical angle of this system could be measured in real experiments. 



PACS numbers: 03.65.Vf, 78.20.Ls 

I. INTRODUCTION 

The low energy effective theory describing three dimen- 
sional topological insulators (TIs) is given by Maxwell 
electromagnetism augmented by an axion like E-B terrrP 
leading to modified constitutive relations. For the bulk 
material this effective theory is valid at energy scales 
below the gap. In the presence of interfaces the mass- 
less modes on the surface would have to be included in 
the effective description of the material and the low en- 
ergy effective description in terms of Maxwell theory with 
modified constitutive relations only applies if the surface 
modes are gapped by an external (time reversal break- 
ing) perturbation. Furthermore the effective description 
is only valid at energy scales below this induced surface 
gap. Such an external breaking can easily be set up ex- 
perimentally, e.g. by a magnetic field. Several poten- 
tial experimental consequences follow from this effective 
theory, such as, for an example, a non-trivial Faraday 
and Kerr rotatiorffl By scanning the external field, the 
topological contribution to the Faraday effect from the 
time reversal breaking field can be cleanly separated from 
the topological contribution. One of the most spectacu- 
lar predictions of this effective theory is the appearance 
of a magnetic monopole mirror charge when solving for 
the static electromagnetic fields sourced by a single point 
charge (located inside a topological trivial insulator such 
as e.g. vacuum) in the presence of a TI interface^. 

For an infinitely extended planar interface, the corre- 
sponding magnetic mirror charge is a pointlike magnetic 
monopole (also carrying some electric charge). As always 
with mirror charges, this monopole of course is not phys- 
ical but simply a mathematical tool to calculate the mag- 
netic fields in the physical region (that is for calculating 
the fields above the interface, the monopole appears to be 
located below the interface and vice versa). Microscopi- 



cally it is surface currents on the interface that source a 
magnetic field with a 1 /r 2 fall-off where r is the distance 
to the mirror monopole. Nevertheless, the magnetic fields 
generated this way are, in the physical region, indistin- 
guishable from the ones generated by a genuine monopole 
and so share some of its properties. In particular, it 
is well known that the electromagnetic fields generated 
by an electric point charge e and a spatially separated 
monopole of magnetic charge g carry a net angular mo- 
mentum which has several interesting properties: it is in- 
dependent of the distance between charge and monopole, 
pointing in the direction from the charge to the monopole 
and proportional to eg (see e.g. Ref. [3]). This total an- 
gular momentum of the composite dyon formed by the 
charge-monopole pair is the sum of the angular momenta 
of the two point particles and this angular momentum 
stored in their fields. Via the spin-statistics theorem this 
shift in the angular momentum of the dyon is often inter- 
preted as a shift in the statistical angle that determines 
the behavior of the multi-dyon wavefunction under the 
exchange of two dyons. For a genuine charge/monopole 
pair, the Dirac quantization of magnetic charges ensures 
that the resulting angular momentum is an integer mul- 
tiple of h/2. So while the statistics of the dyonic system 
can be changed due to angular momentum stored in the 
fields, the net angular momentum is still properly quan- 
tized and so the overall statistical angle is always an in- 
teger multiple of 7r. The full dyonic system is either a 
fermion or a boson. 

As the field of the mirror monopole is indistinguishable 
from the field of a real monopole, the same calculation 
implies that the electromagnetic fields generated by a 
single electric point charge e above the interface of a TI 
carry a non-trivial angular momentum as well. As the 
electric as well as the magnetic field below the interface 
(that is inside the TI) appear to be sourced by coinci- 
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dent charges located at the location of the actual phys- 
ical charge, the contribution to the angular momentum 
from that region of space vanishes. By symmetry, the 
actual angular momentum in the system obtained from 
integrating over the electric and magnetic fields above the 
interface, sourced by the physical electric charge as well 
as the mirror charge, gives exactly half of the angular mo- 
mentum one would get from a genuine charge/monopole 
pair with the same values of e and g (for details see the 
appendix) . However the magnitude of the induced mir- 
ror charge is proportional to the finestructure constant 
a and furthermore depends continuously on the mate- 
rial properties \x and e. Consequently, it generically does 
not obey Dirac quantization conditions. The resulting 
angular momentum is not quantized. As a result, the 
statistics one associates with these charge carriers based 
on their angular momentum is no longer simply fermionic 
or bosonic. Instead they seem to behave as anyons^with 
a statistical angle given bj0 
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Ps = 0/(2tt) is the electromagnetic polarization. It is 
in a topologically trivial material and 1 /2 inside a TI. In 
order to ensure that our expression for &s can be com- 
pared between different unit systems, we explicitly dis- 
played factors of eo and //o- For the bulk of this work 
we'll work with units where eo — Mo = 1- 

The non-trivial statistical angle was interpreted in 
Ref. [5] as a result of the two dimensional nature of the 
TI surface. It is well known that in two spatial dimen- 
sions, anyonic statistic is allowed. The authors of Ref. 2 
proposed that in the vicinity of a TI surface an electric 
point charge indeed turns into an anyon with statisti- 
cal angle 6s and proposed an explicit experimental setup 
that would allow its measurement. This proposal raises 
one important conceptual puzzle: the angular momen- 
tum in the system and hence the inferred 6s is entirely 
independent of the distance zq between the point charge 
and the interface. While it is reasonable to assume that a 
point charge in close vicinity of a TI surface has anyonic 
character, the calculation of &s via the induced angu- 
lar momentum seems to predict that any point charge 
moving freely in three dimensional space would pick up a 
statistical angle #5 provided there is a planar TI interface 
somewhere in the universe at arbitrary large distance zq , 
which is very counter-intuitive and also seems to indicate 
that the anyon is truly 3+1 dimensional in character, con- 
tradicting the fact that 3+1 dimensional anyons should 
be impossible (the exception recently proposed in Ref. 4 
can readily be understood in terms of a more complicated 
topology of configuration space in this case^) . 

One could expect the zg-independence of the angular 
momentum to be an artifact of the special example of 
an infinite planar interface. After all, L/h is dimension- 
less and so could only depend on zq in the form of the 
ratio Zq/cl, where a is another geometric scale in the prob- 
lem. For the infinite plane, no such other scale is present. 



With this puzzle in mind, we analyze the angular mo- 
mentum associated with an electric point charge in the 
vicinity of a TI interface for two different geometries: a 
TI in the shape of a sphere and a TI with a semi-infinite 
tube like shape inside a perfectly conducting cavity. We 
find that in both cases the angular momentum vanishes 
identically for any charge separated from the surface even 
by an infinitesimal amount. These two examples make us 
suspect that the angular momentum will in fact vanish 
for a charge close to (but not right on top) the surface 
of any finite size TI. In hind-sight, this result is not too 
surprising. After all, the microscopic description of the 
topological insulator is in terms of a system of electrons 
and protons with properly quantized charges, obeying 
the standard rules of quantum mechanics. Any state de- 
scribed by this microscopic system has to have a properly 
quantized angular momentum. As long as the effective 
theory correctly captures the long distance behavior of 
the system, it has to obey the quantization conditions 
obeyed by the microscopic constituents. So the angular 
momentum L for a charge close to any compact TI has to 
be an integer multiple of h/2. As a consequence, L can 
not continuously vary as a function of z$ / a. As we expect 
L — >• for zo/a — > 0, it should have been expected that 
L = is indeed the correct answer for all finite values of 
z Q /a. 

At first this result may indicate that the non-trivial 63 
identified in the planar case does not carry over to any 
compact sample and hence would not be measurable. But 
this is too naive. After all, total angular momentum of 
the one-particle system was only taken as a stand-in for 
the statistical angle of the excitations. A more careful 
analysis should directly analyze the two-particle system 
and study the change in action associated with a non- 
trivial loop in configuration space. Performing this anal- 
ysis we find in the case of the planar interface that 63 as 
inferred from the angular momentum in the one-particle 
system only describes the exchange of two particles in 
the limit that the size of the loop I is much larger than 
the separations z^ and z^ of the charges to the surface 
of the TI. This is consistent with the interpretation of 6s 
as a topological effect. For I ~ z^'^ short distance ef- 
fects become important. In this limit details of the path 
matter. But for I 3> d\ 2 the only effect surviving is the 
topological phase (which in this limit is independent of 
the shape of the loop as it should be). So the statistical 
angle governs large loops in configuration space. How- 
ever, for the realistic case of a compact TI of linear size 
a, we should clearly expect significant finite size effects 
in the case that loops are of order the sample size, I ~ a. 
Indeed this expectation is born out. For generic Zq 1 ^ 2 -' 
and a the change in action associated to taking particle 
2 around particle 1 depends crucially on the path and 
has no relation to the 6s obtained from the planar case. 
However, in the intermediate loop size regime 



we once more are able to show that the phase is topolog- 
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ical (independent of shape) and is given by the flat space 
value 6s ■ While not surprising, this analysis clearly lays 
out that any experimental attempt at measuring 6*5 e.g. 
as proposed in Ref. [H has to be set up as to operate in 
this intermediate loop size regime. 

This note is organized as follows: In the next section 
we'll review the solution for the potentials in the case of 
a spherical TI and construct the solution of the poten- 
tials in the case of a tube-like TI. In section 3 we analyze 
the distance dependence of the statistical angle. We first 
show that in the one particle system angular momentum 
vanishes identically in both examples. We then analyze 
explicitly the two-particle system and show, in the case 
of the sphere, that 9s obtained from the flat space anal- 
ysis does govern the intermediate loop size regime. We 
discuss the relevance of these results in section 4. 




II. ELECTRIC AND MAGNETIC POTENTIALS 
IN TWO COMPACT EXAMPLES 



A. Spherical TI 

First let us analyze the electric and magnetic fields for 
a spherical topological insulator and a point-like electric 
charge outside the sphere. The corresponding potentials 
have been worked out in the supplementary material of 
Ref. |U As shown in Fig. [T] a spherical topological insula- 
tor of a radius a and a magneto-electric polarization P3 
is centered at the origin, and a point-like electric charge 
is at (0, 0, d). t\ and [L\ are the dielectric constant and 
the magnetic permeability outside the sphere, £2 and \i 2 
the corresponding quantities inside the sphere. 
Both inside and outside the sphere, the curl of electric 
and magnetic fields is zero, thus we can find a scalar po- 
tentials in both regions: 



E 

(i) 
M 



(2) 



where i = 1,2 stand for inside and outside region. The 
most general solution for the potentials in eq. ^ can be 
written in terms of Legendre polynomials: 



-^E^(COS0)+E^(^)' + 1 ^(COS( 



=£Si(~)'fl(cos*. 

*M =£A(£)'fl(cOS0) 



(3) 



Solving boundary condition for the interface between 
trivial and topological insulator (that is continuity of the 
perpendicular components of D and B as well as the par- 



FIG. 1: A charge q outside a spherical TI. 



allel components of H and E), one arrives at: 

A _ q a 1 r (ei;-e2iWMl + (' + l)/M2]-(2aP 3 ) 2 K*+l) 1 
^ - £l di + ^(2aP 3 ) 2 i(;+l) + (£ 1 (i+l)+€ 2 ;)[VA'i+('+l)/A'2]J 

r> _ q a' r (e 1 l-e 2 l)[l/fi 1 + (l + l)/ fJ , 2 ]-(2aP 3 ) 2 l(l+l) , } 
°l - ei d«+l K2 a P 3 ) 2 i(m) + («i(*+l)+<=2!)[VMi + ('+l)/M2] f J 

r _ n _o}_ (2aP 3 ) 2 i(2; + l) 

°i - 1 d i + i (2aP 3 ) 2 ((m) + (£i(i+l)+e 2 (WA»i+(m)Ai2] 

r>, - -(2aP 3 ) 2 (i+l)(2;+l) u , 

- q d i + i ( 2a p 3) n(l + l) + (e 1 (l + l)+e 2 l)[l/ii 1 + (l + l)/n2\ VV 

The fields here could be considered to be generated by 
a point image electric charge, magnetic monopole, and a 
line of image electric or magnetic charges^. 



B. Semi-infinite rectangular Tl-tube in conducting 
cavity 

Next we want to consider the case of a semi-infinite 
topological insulator (the TI being at z < 0) tube 
with rectangular cross section inside a conducting wall, 
see Fig. [2j An electric point-like charge is located at 
(x',y',z'). As in the spherical case, we can write down 
general forms for electric and magnetic potential: 



E = ^£°mnSin(fc m x)sin(fc„y)e' )w(2: z "> 
+ £ A mn sin(k m x) sin(fc n y)e _7m " z 



^ { E = £ B mn sm(k m x) sm(k n y)e< mnZ 
$ M = £ C mn sa\(k m x) sin(fc„y)e" 7m " J 
*m = £ D mn sin(fc m a;) sin(fc n y)e 7m " z 



(5) 



with O mn being the coefficients of the Green's function 
of a point-like charge in this system: 



sin(fc m x') sm(k„y'). 



(6) 
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Z q(x',v',z') 




FIG. 2: An electric charge q inside the positive half-tube 
z > with dielectric constant ei and magnetic permeability 
/ii in the vicinity of a TI with £2 and ^12 filling the negative 
half-tube z < 0. The tube has a width a, and length b. We 
assume that the walls of the tube are formed by a perfect 
conductor. 



Solving boundary condition, we arrive at: 

2gO mn e-^ z ' (2aP 3 ) 
mn " ( ei +e 2 )(l//i 1 + l//i2) + (2aP 3 ) 2 

Cm 71 — D mn 

r> _ 2gO mra e-^" z '(l/ / i 1 + l/ / i 2 ) 

(ei+e 2 )(l//ii + l/M 2 ) + (2«F 3 ) 2 [> 

f _I , 2(l/^i + l/^) 

1 ex (ei+e 2 )(l/ A11 + l/ At2 ) + (2aP 3 ) 2J 

In this case the geometry of the TI is not really com- 
pact. However, due to the Dirichlet boundary conditions 
on the conducting walls from the Id point of view all our 
gauge fields are massive and exponentially decay at large 
\z\. For all practical purposes, the TI is finite in extend 
in the z-direction. 



III. DISTANCE DEPENDENCE OF 
STATISTICAL ANGLE 

A. Total angular momentum as a global probe of 
the system 

The total angular momentum is given by: 

The extra factor of 1/ (4tt) compared to Ref. |3] is due to 
the fact that we follow Gaussian units as in Ref. [SJ In 
a normal system without a topological insulator, there 
would not be any non-zero angular-momentum of the 
system in the presence of only a static point-like elec- 
tric charge. In the presence of a topological insulator, L 



can be non-zero. In order to address the question about 
distance dependence of the angular momentum and the 
statistical angle in the case of TIs with a compact geom- 
etry, we want to calculate the angular momentum in the 
two examples described in the previous section. 

Let us start with the case of the rectangular tube. 
Plugging our answers from eqs. (J2J), ^ and Q into the 
expression for the angular momentum eq. ^ , we see that 
the integral over x and y can easily be done analytically. 
For example, starting with the x momentum density 

h = y{d x ^Ed y (j)M ~ d y $ E d x <f)M) 
-z(d z $ e d x (j)M - d x (j) E dy(j)M) (9) 

we see that integrating over y first gives an expression 
proportional to 

Jq dyl x ~ mi cos(mi7ra:/a) sin(r7i7rx/a) 

+m cos(m7ra:/a) sin(mi7ra;/a) (10) 

for integers m, mi still to be summed over (as well as 
n and n\ appearing in the coefficients). Further integra- 
tion of this expression over dx then vanishes identically. 
We have confirmed that L y and L z similarly vanish after 
doing the x and y integrals using Mathematica. 

The calculation in the spherical case is a little more 
cumbersome. Plugging eqs. ((2|, ^ and Q into eq. (|8j), 
we can use the recurrence relation of Legendre polyno- 
mials to simplify the integrand. After performing the 9 
integral, the double sums occurring in E x B collapse 
into a single one. Performing the integral in spherical 
coordinates, we once more find that the angular momen- 
tum vanishes identically as long as the electric charge is 
separated even an infinitesimal amount from the surface 
of the sphere, irrespective of the distance d. If we set the 
distance d to zero to begin with and then perform the 
integrals and sums, we get back to the result of eq. 
valid for the infinite half-plane. The fact that for d = 
the sums give back this non-trivial result is a non-trivial 
check. For d/a = one should clearly recover the result 
of the plane, which can also be thought of as a — >• 00. 
What is surprising is that this limit is not smooth. The 
angular momentum vanishes for any finite d. 
As discussed in the introduction, this result should have 
been expected based on microscopic considerations. For 
a real finite system made of electrons and protons total 
angular momentum has to be quantized. So it can not 
continuously depend on d/a. For it to vanish at infinity it 
has to be zero for all d/a. Even for the topological insu- 
lator, genuine fractional angular momentum should only 
exist in the infinite system. With the electron charge be- 
ing on the surface, one can get back the infinite system 
result. 



B. Finite size loop path in the two particle system 

While the angular momentum for a genuine 
charge/monopole pair correctly captures the statis- 
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tical angle, our zero result from the previous subsection 
strongly suggests that in the case of a spatially compact 
TI one should be more careful. To directly obtain the 
statistical angle, the natural thing to do is to look at 
two particles moving adiabatically around each other 
and study the phase change of the action. The action 
of two point particles at a fixed distance Zq above a 
Tl/insulator interface located at z — contains the 
standard terms coupling the point particles to the gauge 
fields: 



S = 



dx^ 
M dr 



(11) 



where the super/sub-scripts (1) and (2) refer to the two 
particles, being the field sourced by particle i, x^ its 
position. Let us for simplicity look at the case where par- 
ticle 1 is kept at the origin in the x-y plane and particle 
2 is taken around a non-trivial loop. 



1. Planar interface 

It is easy to work out the effect of this coupling for 
the planar interface in detail. Note that we can write the 
magnetic field of the mirror monopole of magnetic charge 
g at z = — zq induced by an electric charge e at z = zq 
(and bothat at x = y = 0) in terms of a vector potential: 



4tt 



(1 



COSt 



where (f> and 9 are the angles in a spherical coordinate 
system centered on the location of the mirror monopole 
at (0, 0, — z ). That is the 2nd particle at (x( 2 ),y(2), z o) 
in this coordinate system is located at tan#( 2 ) = P(2)/ z o 



where p 



- y . This is the standard form of the vec- 



tor potential of a monopole with a Dirac string running 
along the negative z-axis. This is the appropriate form 
to use for a mirror charge located below the interface, as 
this mirror charge is only supposed to be used when cal- 
culating fields above the interface (for the magnetic field 
below the interface one would similarly use the A asso- 
ciated to a monopole with a Dirac string running along 
the positive z-axis for a monopole located at z = +zq. 
As there is no physical charge located below the interface 



and so no contribution to eq. ( 11 ) we will not need A in 
this region). 

If we keep one particle fixed at (x — 0, y — 0, z — z ) 
and take the other particle around a loop at a fixed zq 
(4>(2) goes to (j)(2) +2-7T at fixed #(2)), e lS picks up a phase: 



e iS g e JdT<t> {2) A { 4 



1(1- 



)) 



We see that due to the 9(2) dependence already in the 
case of the planar interface, the resulting statistical angle 
depends on path. As tan #(2) = p(2)/ z o, we get different 
answers depending on the size of the loop. When the 
charge is a finite distance zq above the interface and taken 



around a very small loop of size I <C zq we have 9(2) — 
along the whole path and correspondingly do not pick 
up any phase. With a loop large compared to z (that is 
I S> zq), 9(2) 7r /2- In this limit we get a phase e le9 ^ 2 in 
the action independent of the detailed shape of the path 
in the x-y planed This change in the action corresponds 
to e 2tSs where 9g is the statistical angle describing the 
exchange of two particles. Exchanging particle 1 and 2 
twice should correspond to a closed loop in configuration 
space as the one we have been analyzing. So 9s = eg /4. 
Plugging in the g obtained for the mirror charge in Ref. 2 
(rederived in SI units in the appendix here) we do get 
back to eq. 0. 

In conclusion, we find that for a planar interface in the 
realistic case of a finite zq (so that our effective theory 
applies), only a large loop (compared to zq) gives a shape 
independent answer governed by a topological statistical 
angle. For these large loops the statistical angle is half 
of that of a dyon made of a real electron/monopole pair 
in complete agreement with the angular momentum 
calculation in previous section. The modification for 
small loops presumably should be understood as a result 
of short range interactions. 



2. Spherical TI: 

For the case of a spherical TI of radius a we can also 
obtain the vector potential associated with a single point 
charge at (0, 0, d) with d = a + zq. It again only has a <j> 
component with 



A^ — — Coacos( 



1 



sm(9)Pl(cos9). 



The Ci are given in eq. Q. Most importantly, Cq — 
(this would be a net magnetic charge, so it has to van- 
ish). 

We use this expression for A^ and study the phase change 
in spherical coordinates when two particles located at the 
same r = d = a + zq but different 9 are exchanged. We 
take one particle fixed at (0, 0, d) while the other one goes 
from <p to 0+27T at a given 9. The action picks up a phase 
2ttA^. We plot this value for different Zq together with 
the universal 9s in Fig. [3] 

There are a couple of interesting features about this plot. 
First, A^ always vanishes for 9 — > and 9 — > tt, as ex- 
pected after considering the planar case. 
Second, for intermediate 9 between and 7r/2, there is 
a certain non-zero plateau region. For small zq = d — a 
(small compared to a), this plateau region repeats the 9 
behavior of planar case before finite size effect kicks in. 
In this limit the plateau value of A^, goes back to planar 
case result, the universal 9s- 

This can also be seen analytically from the large a asymp- 
totic behavior of A^. To analyze the large radius limit 
it is convenient to rewrite the Ci in a way that makes 
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it explicit that is being sourced by a point mirror 
charge together with a mirror line charge extending from 
the mirror point charge to the origin: 



Ci 



<J2 



Cl 



(-2 



-h + l + i 



-t 2 + l + i 



(12) 



Explicit formulas for point charge gi and the parameters 
t\ 2 an d C12 characterizing the line charge appear in the 
supplemental material of Ref. H] for e = /i = 1. Most 
importantly, the point charge (72 goes to the planar value 
g in the d — > a limit. We have confirmed with Mathe- 
matica that the same decomposition holds for arbitrary 
e and \i. The general formulas for 32, ti,a an d Ci t 2 are 
too lengthy to reproduce here. Taking a large at fixed 
d, one sees that every term in the potential scales as 1/a 
and so the potential seems to vanish in the large a limit. 
This conclusion is too fast as the infinite sum can alter 
the behavior. At = (where P;(l) = 1 for all I) the 
sum can be performed analytically and we see that the 
r = a contribution for the point charge sums up to the ex- 
pected gil(d— a) = g2/zo, that is it remains finite in the 
large a limit. As the sum over Legendre polynomials by 
construction just represents a standard \/r Coulomb po- 
tential, the point charge term automatically reproduces 
the planar contribution to 63- At 9 = the contribution 
of the line charge can also be summed up. The resulting 
Hypergeometric function vanishes as 1/a in the large a 
limit. As Pi < 1 for 9 7^ it is clear that at non-zero 
9 the line charge contribution has to vanish at least as 
fast as for 9 = 0. So in the a — > 00 limit, the line charge 
contribution can be neglected compared to the point-like 
one. This result has also been confirmed numerically in 
Ref. H 

This shows that for loops of size I <C a we recover the 
planar result. However in the planar case we found that 
for I 3> zq the exchange is governed by the universal 9$. 
From this we conclude that for a compact geometry it is 
the intermediate size loops with zq <C I -C a which arc 
governed by the universal topological phase of eq. ([lj. 
While we only explicitly demonstrated this result in the 
case of a sphere, we believe it to be true in general. 



IV. DISCUSSION 

We find zero angular momentum for a charge outside 
a spatially compact TI in two examples. While obvious 
from the microscopic point of view of electrons and pro- 
tons, the result looks somewhat surprising starting from 
Maxwell's equations of a topological insulator; in fact we 
haven not been able to give an analytic proof based on 
this effective theory that the total angular momentum 
always vanishes in a finite system even though the mi- 
croscopic point of view strongly suggest that this is true. 
By looking at the two particle system we find that the 
universal statistical phase governs the behavior of inter- 
mediate size loops. In realistic experiments, our results 
indicate that to observe such an anyon in the topological 
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FIG. 3: A$ versus 6 with different zo's for the mirror charge 
on an electron, zq = 0.05a for the dotted line and zo = 0.1a 
for the dashed line and zq = 0.0005a for the solid curved 
line. The horizontal straight line stands for the universal angle 
8s- The thickness of the zq = 0.0005a line represents the 
numerical uncertainty due to the slow convergence of the large 
I truncation. In this plot we chose a — 1 and all permittivities 
and permeabilities equal 1. The small numerical value of 9s 
is due to the appearance of a in the expression for the mirror 
charge. 



insulator set-up as suggested in Ref. 21 one should use a 
large bulk sample with a small superconductor loop close 
to it, to get non-zero flux. The size of the superconductor 
loop has to be in the intermediate regime we identified, 
that is much smaller than the size of the sample but much 
large than the distance between loop and sample. Our 
results are consistent with the 2+1 dimensional nature 
of the anyons, as the universal statistical angle ceases 
to accurately describe the two particle systems once the 
particles are removed from the surface of the TI beyond 
a distance of order the sample size. 
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Appendix A: Derivation of the statistical angle in SI 

units 



Consider a planar interface between a topological in- 
sulator (with non-trivial /12 and 62 as well as 9 = ir, that 
is P3 = 1/2) and a trivial insulator (with ei and Hi). 
Consider a single static point charge q inside the trivial 
material a distance zq away from the surface of the TI. 

We want to find the mirror charges for general fj, and 
e. Let the TI occupy the z < region of space. Writing 



Maxwell's equation as usual as 



V-D = p e ,WxH 



dD 
~dt 



v4 * dB f 



7 



The constitutive relations in SI unites are 



D = eE - e Q a- (c B), 

7T 



c H 



c B 



a- 



e e 

i , 

'?r ho 



We can introduce potentials $ e ,m with E = — V$ e 
and B = — V<& m . Above the interface they are given by 
(note that we weight all electric mirror charges by eo for 
convenience; we also use e\ for q e on both sides of the 
interface; these are just definitions of our mirror charges) 



(2) 



47reii?i 47reoi?2 



(2) 



4ttR 2 



and below by 



= ge + gg 

47reii?i 47re -Ri 

,(1,2) 



o (1) 
47ri?i 



where <7e,m are the mirror charges locate a distance <i 

above (for qi]L) or below (for gifm) the interface. Rf = 
x 2 +y 2 + (z -z) 2 andi? 2 . = x 2 +y 2 + (z + z) 2 . Maxwell's 
equations in the absence of surface currents or charges as 
usual demand continuity of D±, B±, H\\ and E\\. As, at 
z = 0, R\ = J?2, 9 Z -Ri = —d z R,2 this demands: 



- ^ 2) ) 


= ( — 9e " 






™ 


= -g( 2 ) 
ym 


a (2) 


a (1) 






Mi 


M2 
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g e /ei + qi^/eo 

i 

7T /X C0 



From this the mirror charges can easily be found: 



T (2) 



and 



7 (2) 



m = 1 2 ^fg 

/m C o( £l + e2 )(^ + ^) + eo a 2 i 



n) _ e (ei - e a )(ff + ff) - e a^ 



e a 



2&L 



The system consisting a the charge g and a magnetic 
charge g m gives rise to an angular momentum 

f Vim . 

where f is the unit vector pointing from the electric to 
the magnetic charge. Clearly L vanishes when the two 
charges are sitting on top of each other (as in this case E 



and B are parallel, so the momentum density and hence 
the angular momentum density vanish identically). 

For the interface, we need to calculate the contribu- 
tions to the angular momentum in the two regions in- 
dependently. Inside the TI both electric and magnetic 
fields are pointing radially outward from the point at 
z = +zq, so the angular momentum vanishes (again, E 
and B are parallel and so the Poynting vector vanishes 
identically). For the region above the interface, we get 
a non-zero contribution to the angular momentum due 
to the charge/monopole system formed by the original 

charge q at z — +z and the mirror magnetic charge q m ' 

(2) 

at z = — zo (the electric mirror charge q\ at z = — zq 
does not contribute, as the charge at z = +zq is purely 
electric). If we calculate the angular momentum with 
respect to the origin x = y = z = 0, we see that the 
integrand 



rx(ExB) 



rx[(r-z e z )x(f+z e z )] 



\r- 



-z e, 
rx Irxe 



| 2 |?+zo<M 2 



Z \f-z e z \ 2 \f+z e z \ 2 



(Al) 



is symmetric under r — > —f. So we get equal contribu- 
tions to the angular momentum from the lower and the 
upper half plane. As in our case we only get a contribu- 
tion from the upper half plane, the angular momentum to 
the charge/mirror-charge system is exactly half of what 
it would be for a genuine charge/monopole pair, that is 
(setting 9 = ir, q = — e in the expressions above) 



1 ae~ 

8tt - 47TCO ( £l+e2 )(m + g) +eoQ 2 

g£h 

- (£i + £2 )( ™ + ™ )+Q 2 



(A2) 



where in the last step we used the standard SI definition 
for a, 



47re fic 



Note that for a given q, say q = — e, but for general n 
and e this is certainly not going to obey the Dirac quan- 
tization condition. The statistical angle one would want 
to associate with a dyon with angular momentum L z is 
9s = 2nL z /h, so that the statistical angle 9s is for inte- 
ger spins (bosons) and n for half-integer spins (fermions). 
For a charge — e in the presence of a TI surface we there- 
fore obtain 



\6s\ = 2tt- 



TsL + e)(£a + — ) + a 2 ' 
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6 Our expression differs by a factor of 2-k from the one quoted 
irP. A factor of 4-zr seems to be due to the fact that irP, 
which works in Gaussian units, the formula for the sta- 
tistical angle has been normalized to yield 9s — ir for 
a monopole of flux hc/e. However latter has been taken 
to correspond to g = hc/e, which would be correct if 
V • B — g8(r) as it is in SI units. However in Gaussian 
units V • B = 4ivg5(r) and so the formula for 9 s irP in 
terms of general e and g should be modified by a factor of 
4-7T. The remaining factor of 2 mismatch seems to be due 
to the fact that irP the factor of 1/2 from the fact that 



only the half of space above the interface contributes to the 
angular momentum has not been taken into account. For 
the reader's convenience we rederive the expression of the 
statistical angle in SI units in the appendix, where interme- 
diate steps can be readily compared to standard textbook 
expressions from e.gP. This confirms the factor of 2-k as 
we have it here. For the bulk of the paper we follow^ and 
use Gaussian units. The final expression for the statistical 
angle has to be independent of choices of units when only 
expressed in terms of the fine structure constant a as well 
as the ratios e/eo and n/fj,o. 
7 These answers also obvious when one considers that V x 
A — B, so the line integral of A we are doing is equal to 
the total magnetic flux through the surface enclosed by the 
loop. For 0(2) = the loop has vanishing area and no flux. 
For 9(2) = 7r/2 the loop captures all the flux in the northern 
hemisphere. 



